Non-vanishing spin Hall currents in disordered spin-orbit coupling systems 



O 

o 

(N 

> 
O 

(N 



K. Nomura, 1 Jairo Sinova, 2 T. Jungwirth, 3 - 4 ' 1 Q. Niu, 1 and A.H. MacDonald 1 

department of Physics, University of Texas at Austin, Austin TX 78712-1081, USA 
2 Department of Physics, Texas A&M University, College Station, TX 77843-4242, USA 
'^Institute of Physics, ASCR, Cukrovarnickd 10, 162 53 Praha 6, Czech Republic 
* School of Physics and Astronomy, University of Nottingham, University Park, Nottingham, NG7 2RD, UK 

(Dated: February 2, 2008) 

Spin-orbit coupling induced spin Hall currents are generic in metals and doped semiconductors. 
It has recently been argued that the spin Hall conductivity can be dominated by an intrinsic contri- 
bution that follows from Bloch state distortion in the presence of an electric field. Here we report 
on an numerical demonstration of the robustness of this effect in the presence of disorder scattering 
for the case of a two-dimensional electron-gas with Rashba spin-orbit interactions. 

PACS numbers: 72.10.-d, 72.15.Gd, 73.50.Jt 
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Semiconductor spintronics research over the past 
decade has concentrated on the properties of spin- 
polarized carriers created by optical orientation, on 
the search for new ferromagnetic semiconductors with 
more favorable properties, and on the injection of spin- 
polarized carriers into semiconductors from ferromag- 
netic metals Q, H Q • There has recent ly b een a fl urry 
of theoretical interest 0JLE 0, EM M, El 13 111 
Q E1EI El El El IflM Slip in the spin Hall 
effect 313,11, i.e. in transverse spin currents induced 
by an electric field. Murakami et al. and Sinova Q 
et al. have argued in different contexts that the spin 
Hall conductivity can be dominated by a contribution 
that follows from the distortion of Bloch electrons by 
an electric field and therefore approaches an intrinsic 
value in the clean limit. This conclusion has recently 
been questioned, for the case of two-dimensional elec- 
trons with Rashba spin-orbit interactions in particular, 
by several authors El El El El El E3 motivated by 
a number of different considerations, some of which are 
related to controversies [M El El El El El El that 
have long surrounded the theory of the anomalous Hall 
effect in ferromagnetic metals and semiconductors. In 
this Rapid Communication we report on a study based 
on numerically exact evaluation of the linear-response- 
theory Kubo-formula expression for the spin Hall con- 
ductivity. We demonstrate that the intrinsic spin Hall 
effect is robust in the presence of disorder, falling to zero 
only when the life-time broadening energy is larger than 
the spin-orbit splitting of the bands. The correlations be- 
tween spin-orientation and velocity in the presence of an 
electric field that lie behind the intrinsic spin Hall effect 
are not diminished by weak disorder. 

We consider a two-dimensional electron system with 
the Rashba spin-orbit interaction(R2DES): 



H = p 2 /2m + A[p x z] • er/H 



V. 



(1) 



where er is the Pauli matrix, m is the effective mass, 
and A is the Rashba spin-orbit coupling constant. When 
the disorder potential V in Eq. (JTJ is absent, p = Kk 
is a good quantum number. The Rashba spin-orbit in- 
teraction term can be viewed as Zeeman coupling to a 



k-dependent effective magnetic field A = (2A)z x k. 
The V = eigenstates are therefore the S = 1/2 
spinors oriented parallel and antiparallel to these fields: 
|k±) = [ =Fie - ^, 1 ] e l ^' r /\/2fT, and the two eigenvalues 
at a given k are split by 2A|k|. Here <j> = tan -1 (k x /k y ), f2 
is the system area and we have applied periodic boundary 
conditions. As explained in Ref.yJ, an electric field in the 
x-direction causes Rashba spinors to tilt out of the x-y 
plane giving rise to an intrinsic spin Hall effect. The key 
issue in dispute is whether or not the velocity-dependent 
spinor tilts vanish when quasiparticle disorder scattering 
is properly taken into account. To address this subtle 
issue without making any assumptions which might prej- 
udice the conclusion, we evaluate the Kubo formula for 
the spin Hall conductivity using the exact single-particle 
eigenstates of a disordered finite area two-dimensional 
electron system with Rashba spin-orbit interactions. 

Our disorder potential consists of randomly centered 
scatterers that have strength uq and a Gaussian spa- 
tial profile with range l v . The potential matrix elements 
satisfy Kk£7|V|kV)| 2 = [mul/tyS^, exp(-|k - k'| 2 Z 2 ), 
where the density of scatterers nt (intended to represent 
remote ionized donors) is set equal to the electron den- 
sity. It is widely recognized that 2DES disorder poten- 
tials can have long correlation lengths up to ~ 100 [nm]. 
To examine how our conclusions depend on the range of 
the disorder potential, we have performed calculations for 
correlation lengths ranging from l v ~ to l v ~ 100[nm]. 

We diagonalize the Hamiltonian in the A = eigen- 
state representation and introduce a hard cutoff at a suf- 
ficiently large momentum A. For a fixed particle density, 
the number of electron N e and the system size are re- 
lated by fl = L 2 = N e /n e . Our conclusions are based 
on calculations with N e up to 2258. For n e = 0.6 x 10 11 
[cm -2 ] the system size is up to L = 2 [/im], longer than 
the characteristic microscopic length scales, the mean- 
free path (I - 10 2 - 10 3 
(X F = 2n/k F = 101 
range (l v < 100 [nm]). The system size in these simula- 
tions is comparable to that of typical 2DES channels in 
electronic devices. We fix the effective mass at the bulk 



nm]), the Fermi wavelength 
nm]), and the disorder potential 
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GaAs value, m = 0.067m e , where m e is the bare electron 
mass and perform calculations over a wide range of A and 
uq values. 

The Kubo formula expression for the z spin component 
of the spin Hall conductivity is: 



1 
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labels ex- 
and the charge 
dH/dp 



where f(E) is the Fermi function, n 
act eigenstates with eigenvalues E n 
and spin current operators are j = — e 
-e(p/m + Az x tr/h) and j 2 = {dH/dp, %a,}/2 = 
pa z /m respectively 8] . In finite size calculations the 
electric field turn on time r\~ x must be shorter than the 
transit time in the simulation cell in order to obtain the 
correct thermodynamic limit for the conductivity. In the 
metallic limit of interest here, rj must exceed the simu- 
lation cell level spacing but be smaller than all intensive 
energy scales. In the dc u = limit, is real with a 
dissipative contribution that comes from the ir\ term in 
the denominator and a reactive contribution that comes 
from the imaginary part of the matrix element product. 

Typical numerical results for the disorder and spin- 
orbit coupling strength dependence of the spin Hall con- 
ductivity OsH — o% y (uj — 0) are illustrated in Fig.l. 
(These calculations are for l v ~ 80 [nm].) We find that 
in the strong Rashba coupling, weak-disorder regime the 
spin Hall conductivity is close to the (universal) intrin- 
sic value for this model, and that it decreases for weaker 
spin-orbit coupling and stronger disorder. Experimen- 
tally, Rashba spin-orbit coupling strength can be varied 
over a wide range by tuning a gate field |33Ll34j . We have 
varied the spin-orbit coupling strength at the Fermi en- 
ergy Xk F from 0.1e F to 0Ae F . The system size for the cal- 
culations summarized by Fig.l was 1500nm. The range 
we have chosen for disorder strength values was based on 
the golden-rule expression for the transport scattering 

rateHl, H/t = 2tt^' \V(k-k')\ 2 (l-k-k )5{e k > -e F ). 
The golden-rule combined with Boltzmann transport the- 
ory yields the Drude expression for the longitudinal con- 
ductivity, <td = ne 2 r/m = 2e F T(e 2 /h). Using these ap- 
proximate estimates, we have varied the disorder strength 
so that €pT covers the range 2 — 20, typical for two- 
dimensional electron systems. For GaAs materials pa- 
rameters, the disorder strength range that we consider 
corresponds to mean-free paths Z = 70 — 700 [nm] . We 
note that in the case of short-range scatterers (l v ~ 10 
[nm]) the transport lifetime r defined above is not so dif- 
ferent from the momentum lifetime To given by H/tq = 
27r Ek' l^( k " k ')| 2 ^(e k ' - e F ) (l v ~ 10 [nm]), whereas 
these quantities differ substantially for longer (and more 
realistic) correlation lengths. In what follows we take 
h = 1 so that r" 1 has energy units. These results demon- 
strate that for this model <j s h is to reasonable accuracy 
a function of only Xk F r, the ratio of the spin-orbit split- 
ting to the quasiparticle state lifetime broadening. The 
intrinsic spin Hall conductivity survives provided that 
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FIG. 1: Spin Hall conductivity a gH as a function of cft 
and \k F /e F at e F = 2.15 [meV] and n e = 0.6 x 10 11 [cm -2 ]. 
For these calculations the system size is L = 1500 [nm] and 
l v = 80 [nm]. Note that the conductivity depends mainly on 
XkpT and that, because our interest is limited to the metal- 
lic regime, our calculation range does not address the strong 
scattering limit r — > 0. 



Xk F T > 1. 

Fig. 2 illustrates some typical system size dependences 
of the finite-size longitudinal a xx and spin Hall <j s h con- 
ductivities. The size-dependence of transport coefficients 
in disordered systems can reflect quantum corrections 
to Boltzmann transport theory due to the interference 




FIG. 2: Left: Size dependence of the longitudinal conduc- 
tivity a xx as a function of L/l for A = 0, e F r — 2.0, and 
l v = 20 [nm]; Middle and Right: L/l dependence of the spin 
Hall conductivity asH for l v = 20 [nm]) and Xk F /e F — 0.3. 
The middle panel is for a strongly disordered system in which 
e F r = 1 while the right panel is for the a weakly disordered 
system in which e F r — 5. 
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effects that cause localization. In two-dimensions, scal- 
ing theory and microscopic perturbative calculations pre- 
dict a xx corrections that depend on spin-orbit coupling 
strength and can grow when the system size L is larger 
than the mean- free path I. The conductivity is expected 
to decay exponentially with system size in the strongly lo- 
calized region. Numerical a xx results for the strongly 
disordered case epr = 2, A = 0, and l v — 20 [nm], shown 
in the left panel of Fig. 2, are consistent with expectations 
for this thoroughly studied quantity. |35j ] Our main inter- 
est at present, however, is the system size dependence of 
the spin Hall conductivity a s u and particularly in estab- 
lishing whether or not it vanishes in the limit L — ► 00. 
For <j s h , L should be compared with both I and with the 
spin-orbit length L so — 1/ (Xkpr). In the middle panel of 
Fig. [[2] L so w 31 is the longer intensive length scale, with 
some system size apparent up to L/L so ~ 10. For the 
more weakly disordered case in the right panel I is longer 
and no systematic L/l dependence was found. These nu- 
merical results appear to establish rather unambiguously 
that liniL^oo a sH ^ 0. 

The intrinsic spin Hall effect in the R2DEG is due to 
a correlation [3] between quasiparticle velocity and the 
z-component of spin induced by an electric field; for an 
electric field in the x-direction, an up spin is induced in 
positive y-component velocity majority-band states and 
a corresponding down spin at negative velocities. After 
summing over bands, coherence is confined in momentum 
space to the annulus of singly-occupied states. These re- 
sponses are induced by the interband matrix elements of 
the perturbation term in the Hamiltonian that accounts 
for the spatially uniform electric field. Since the observ- 
able we are interested in here, the spin Hall current, is 
purely off-diagonal in band indices, its response depends 
on interband coherence alone and not at all on the altered 
Bloch state occupation probabilities that dominate most 
transport coefficients in metals and are the focus of Boltz- 
mann transport theory. If the spin Hall conductivity were 
to vanish because of disorder scattering, the intrinsic in- 
terband coherence would either have to be cancelled at 
all wavevectors, or be cancelled by stronger coherences 
induced in a narrow transport window (presumably of 
width 1/t) centered on the Fermi circles. 

In Fig. 3 we compare the exact linear-response 
momentum-dependent ^-direction spin-density (and 
hence interband coherence) for a disorder-free system 
(left panel) with Xkp/ep — 0.2 with that of a disor- 
dered system (right panel) with the same spin-orbit in- 
teraction strength and epr = 3.2. ( l v /Xp = 0.2 for 
the calculations illustrated in Fig. 3.) Both quantities are 
proportional to the electric field and are plotted in the 
same units. These results were obtained from the same 
linear response theory expressions used in Eq.(2) with 
S,(k) = */2 |ka)(ka| = (|k+)(k - | + |k-)(k+ |)/2 
substituted for the spin current j*. The disorder aver- 
aged spin Hall conductivity and longitudinal conductivity 
in this case are a S H/(e/8n) — 0.64 and <J xx /(e 2 /h) = 5.1 
at epT — 3.2. Our numerical calculations demonstrate 
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FIG. 3: Electric field induced spin distribution S z (k)/eE as 
a function of wave vector for the clean limit (left side) and for 
an eft = 3.2 I v Xf = 0.2 disorder model (right side). 



that the coherence is not changed qualitatively by impu- 
rity scattering, maintaining the same angle dependence 
as it is spread in momentum space. In particularly there 
is no evidence that the direction averaged coherence is 
either cancelled uniformly or cancelled by a strong con- 
tribution more narrowly centered on the two Fermi cir- 
cles. 

The subtleties that confuse theories of the spin Hall 
conductivity in a R2DES are related to issues that arise 
quite generally in the linear-response theory analysis of 
non-dissipative tr ansp ort coefficients, like the anomalous 
Hall conductivity [37| of a ferromagnet, the ordinary Hall 
conductivity of a paramagnet, and the spin Hall con- 
ductivity of other paramagnetic metals. From an exact 
eigenstate Kubo formula point of view, these transport 
coefficients can be dominated by reactive contributions 
that come from states far from the Fermi level and are not 
associated with electric field induced level crossings and 
dissipation. In the spin and anomalous Hall effect cases, 
the reactive contributions do not vanish in the limit of 
a perfect crystal, instead approaching an intrinsic value. 
The currents accounted for by these intrinsic Hall coef- 
ficients can be viewed as corresponding to equilibrium 
currents that flow in an effective periodic systems whose 
symmetry has been reduced by the electric field. This 
point has been emphasized recently by Rashba 23] , who 
argues on this basis that the intrinsic response is a tran- 
sient that will be attenuated within a relaxation time r 
scale after the electric field is turned on. Similar argu- 
ments have been made concerning the intrinsic contri- 
bution to the anomalous Hall effect. [2^] The specific in- 
stance studied here is perhaps an especially simple exam- 
ple of this class of effects, precisely because S z (k) and the 
spin Hall current are purely off-diagonal in band indices. 
We conjecture, as an extrapolation from the present nu- 
merical study, that the part of the density-matrix linear 
response that is off-diagonal in band index always ap- 
proaches its intrinsic value in the weak disorder limit. 
The spin Hall current operator, like the charge current 
operator in the case of the anomalous Hall effect, will 
also have intraband matrix elements in the general case. 
We expect that these can in general lead to extrinsic in- 
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traband contributions to the linear response conductivity 
that remain finite in the weak disorder scattering limit. 

In a realistic sample with boundaries, spin density is 
accumulated at the sample edge by the spin currents. We 
expect that edge spin accumulations can be measured ex- 
perimentally. Stevens et al^j^ have recently reported on 
a remarkable optical measurement of accumulation due 
to non-linear response spin currents using a spatially re- 
solved pump-probe technique in GaAs/AlGaAs quantum 
wells. Similar luminescence polarization measurements 
should be able to detect electrically generated linear re- 
sponse spin Hall currents. 

In summary we calculated the spin Hall conductivity 
in a disordered system with Rashba spin-orbit coupling 
using the exactly evaluated eigenstates of the Hamilto- 
nian and the Kubo linear response theory. We find that 
the field induced spin Hall current of this model ap- 
proaches its intrinsic value in the limit of weak disorder 



scattering. 
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Note added. — After this work was completed and sub- 
mitted several preprints appeared reporting on related 
numerical simulations [3^, |3j| °f s P m Hall conduc- 
tance in finite samples with contacts. These studies reach 
similar conclusions on the robustness of spin Hall effects. 
Very recently two experimental preprints |4lll42j | have ap- 
peared which report detection of edge spin accumulation 
due to spin Hall currents. 
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